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The transport properties of diluted magnetic semiconductors (DMS) are calculated using dynam- 
ical mean field theory (DMFT) and Boltzmann transport theory. Within DMFT we study the 
density of states and the dc-resistivity, which are strongly parameter dependent such as tempera- 
ture, doping, density of the carriers, and the strength of the carrier-local impurity spin exchange 
coupling. Characteristic qualitative features are found distinguishing weak, intermediate, and strong 
carrier-spin coupling and allowing quantitative determination of important parameters defining the 
underlying ferromagnetic mechanism. We find that spin-disorder scattering, formation of bound 
state, and the population of the minority spin band are all operational in DMFT in different param- 
eter range. We also develop a complementary Boltzmann transport theory for scattering by screened 
ionized impurities. The difference in the screening properties between paramagnetic (T > T c ) and 
ferromagnetic (T < T c ) states gives rise to the temperature dependence (increase or decrease) of 
resistivity, depending on the carrier density, as the system goes from the paramagnetic phase to 
the ferromagnetic phase. The metallic behavior below T c for optimally doped DMS samples can 
be explained in the Boltzmann theory by temperature dependent screening and thermal change of 
carrier spin polarization. 

PACS numbers: 75.50.Pp, 75.10.-b, 75.30. Hx, 75.20.Hr 



I. INTRODUCTION 

Diluted magnetic semiconductors (DMS)i2i2i4 with 
transition temperatures as high as 150K (GaMnAs at 
5% MniiSii) or even above room temperature (GaMnN, 
GaMnP— ) are attracting much attention lately in part 
because of possible 'spintronic' applications^. The proto- 
typical DMS material is Gai-^Mn^As with the Mn ions 
substitutionally replacing Ga at the cation sites. It is 
widely believed that the ferromagnetism in this mate- 
rial is carrier induced with holes donated by Mn ions 
mediating a ferromagnetic interaction between the ran- 
domly localized Mn spins&. The coupling of carrier spins 
(holes in GaMnAs) and localized moments (Mn impuri- 
ties) gives rise to the unique magnetic and transport (as 
well as optical) properties in DMS. 

DMS transport properties are influenced by the ex- 
change interaction between the carriers and the localized 
moments as spin fluctuation scattering contributes to the 
resistivity. The experimentally measured dc-resistivity in 
DMS materials^iSi shows interesting behavior strongly 
depending on the concentration of the magnetic impurity 
and temperature. In Im-ajMn^As^ and Cai^Mn^A*^ 
with low Mn concentration (x < 0.03) only an insulat- 
ing behavior has been observed in transport measure- 
ments. However, near optimal doping x ~ 0.05, where 
the highest value of T c is reported, the non-monotonic 
behavior (insulator-metal-insulator) of the resistivity as 
a function of temperature is observed. A resistivity peak 
appears near the critical temperature (T c ), and the re- 
sistivity shows metallic behavior (dpdc/dT > 0) below T c 
and insulating behavior (dpdc/dT < 0) at higher tem- 
peratures. The peak has been understood as the critical 
scattering effects at T c of spin fluctuations^. 

In this paper we present theoretical calculations for 



DMS transport properties and study the role of the 
carrier-spin coupling which are crucial for ferromagnetic 
properties. Transport measurements have proven useful 
in understanding the physics of the colossal magnetore- 
sistance (CMR) manganites where carrier-spin coupling 
is also crucial^. In our calculation we use a recently de- 
veloped non-perturbative method, the "dynamical mean- 
field theory" (DMFT)jii for calculating the DMS trans- 
port properties. A non-perturbative method is needed 
in DMS materials because the crucial physics involves 
bound-state formation and other aspects of intermedi- 
ate carrier-spin couplings not accessible to perturbative 
methods. (We note that the most interesting phenomena 
in DMS involve intermediate couplings and intermediate 
temperatures. This regime is very difficult to treat by 
standard analytical or numerical methods.) The DMFT 
has been recently applied to the DMS system to calcu- 
late the magnetic transition temperature and the optical 
conductivitjiiSii^. DMFT is essentially a lattice quan- 
tum version of the Weiss mean field theory where the 
appropriate density of states (including impurity band 
formation) along with temporal fluctuations are incorpo- 
rated within an effective local field theory. An important 
ingredient of DMFT— is that it reduces to the standard 
RKKY physics in the weak-coupling regime and the dou- 
ble exchange physics in the strong coupling regime. 

Our DMFT results show interesting dependence of re- 
sistivity on carrier-spin coupling (J), carrier density (n), 
doping of the magnetic material (x), and temperature 
(T) revealing key features of the underlying physics. We 
find that our results show many similarities to colossal 
magnetoresistance (CMR) manganitesiS, especially, for 
a large carrier-spin coupling and near half filling of the 
impurity band (i.e. in the double exchange regime) since 
the spin-disorder scattering dominates in this parameter 
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range. However, the DMS dc-resistivity exhibits novel 
features not found in the CMR in the weak coupling 
RKKY limit (J < 1.0*, where * is the band width of 
the carrier). These novel features arise mainly from the 
bound state formation of the carrier-local spins and car- 
rier occupation of the minority spin band. The formation 
of the bound state gives rise to an insulating behavior 
even if the carrier localization effects are not taken into 
account. Experimental observation of our predictions 
should lead to crucial information about bound state for- 
mation and impurity band physics in this problem. Even 
though we focus on III-V compound based DMS such 
as Gai-^Mn^As, the results presented in this paper are 
general for all DMS materials. 

An important limitation of DMFT is that it is a non- 
perturbative local theory that cannot really incorporate 
spatial fluctuations, and therefore resistive charged im- 
purity scattering with its strong momentum dependence 
is essentially impossible to handle in DMFT. We there- 
fore consider the Boltzmann transport theory to calcu- 
late the resistivity of the DMS systemsi^. The Boltz- 
mann theory is used to calculate the charged impurity 
disorder limited DMS resistivity (whereas the DMFT is 
used for obtaining the spin disorder limited DMS resis- 
tivity) . Charged impurity disorder arises here both from 
the ionized Mn acceptor in GaMnAs as well as other 
charged defects/impurities invariably present in a semi- 
conductor. Using relaxation time approximation we as- 
sume that the Boltzmann resistivity is due to ionized 
impurities. The carriers (holes) are scattered by the 
screened Coulomb potential, which we calculate using 
the linearized Thomas-Fermi (TF) approximation and 
the random phase approximation (RPA). We find that 
the dominant temperature dependence of the resistivity 
comes from the change in the screening length in the 
high density metallic samples. In the metallic GaM- 
nAs DMS samples the change in the measured resistiv- 
ity, when the temperature goes from T c to zero, is about 
20% in good agreement with our calculation. In the fer- 
romagnetic state (T < T c ), as the temperature decreases 
from T — T c , we find strong temperature dependence of 
the resistivity arising from the low temperature screening 
function and the thermal change of the carrier densities in 
each spin-split subband. When all the carriers are polar- 
ized (this happens in the low density limit) the screening 
function is almost independent of temperature and the 
resistivity is also temperature independent. Inclusion of 
both spin disorder and ionized impurity disorder in DMS 
transport is the important ingredient of our theory. 

This paper is organized as follows. In Sec. II, we de- 
scribe our model and our theoretical approach based on 
the dynamical mean-field theory (DMFT). In Sec. Ill, 
we study, in detail, the density of states and describe the 
formation of the spin polarized impurity band. In Sec. 
IV, we provide the results of calculated DMS resistivity 
within DMFT. In Sec. V we calculate the transport resis- 
tivity within the Boltzmann transport theory for ionized 
impurities. In Sec. VI, we summarize our qualitative 



findings and providing a critical discussion of the appli- 
cability of our results to DMS systems. A brief conclusion 
is given in Sec. VII. 



II. MODEL AND FORMALISM 

Our basic model of DMS systems is that of magnetic 
dopants ("impurities") interacting through a local ex- 
change coupling with carriers in the host semiconductor 
material. The generally accepted Hamiltonian of the sys- 
tem is given by 



H = H, 



host 



H 



M 



H AF , 



(1) 



where -Hh os t describes carrier propagation in the host 
semiconductor band. For simplicity we consider a host 
material with a single non-degenerate band. We there- 
fore write 



host 



E 



2m 



+ V R (x) 



1>a(x), (2) 



where a is the spin index and Vr is a random potential 
arising from non-magnetic defects in the material (e.g., 
As antisite defects, unintentional background charged 
impurities, etc.). The second (magnetic) term in Eq. 
||TJ, Hmi describes coupling of the carriers to an array of 
impurity (e.g. Mn) spins at positions i?;, 
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where J is the local exchange coupling between the spin 
of the magnetic impurity and the the spins of the semi- 
conductor carriers, W is the (Coulombic) potential aris- 
ing from the magnetic dopant, Ri are the positions of 
the magnetic dopants, and a is the Pauli matrix. Here 
we absorb the magnitude of the impurity spin into the 
coupling J (which we take to be positive), and repre- 
sent the spin direction by the unit vector S. The third 
term in Eq. ifTJl. Haf, is the direct Mn— Mn short-range 
antiferromagnetic exchange interaction 



Haf — £ Jaf{Ri — Rj)S{Sj 



(4) 



1,3 



where Jaf is a direct antiferromagnetic exchange cou- 
pling between impurity spins. 

In this section we approximate our model by neglect- 
ing the nonmagnetic random potential, Vr, and the di- 
rect antiferromagnetic exchange interaction, Jaf ■ Lat- 
tice defects may be playing an important role in deter- 
mining magnetic and transport properties of the samples, 
but we assume here that these defects enter our theory 
only in determining the basic parameters of the model, 
namely, the density of magnetically active dopants rij, the 
hole density n c , and perhaps the local effective exchange 
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coupling J between the holes and the magnetic impuri- 
ties, and do not include any defects into our model ex- 
plicitly. We include charged impurity scattering through 
the Boltzmann equation in section IV of this paper. We 
believe that the effects of the antiferromagnetic cou- 
pling between magnetic impurities are either negligibly 
small or incorporated into the effective parameters of the 
model. Actually, in the parameter range of interest to 
us (x <C 1), where DMS ferromagnetism typically oc- 
curs, the magnetic impurities are separated from each 
other by non-magnetic atoms, and this short-range an- 
tiferromagnetic interaction, which rapidly decays with 
the distance, should be negligible. These approximations 
are nonessential and are done in the spirit of identifying 
the minimal DMS magnetic model of interest. Both of 
these effects, which may be of quantitative importance in 
some situations, can be included in the theory by adjust- 
ing the parameters of the model or perhaps at the cost 
of introducing more unknown parameters characterizing 
these interactions. Recently several theories for DMS fer- 
romagnetism explicitly including spatial disorder effects 
and antiferromagnetic coupling have been developed^. 

In our DMFT DMS model there are two sources of 
coupling between the carrier and the impurity magnetic 
moment: a spin-spin coupling (J) and a potential scatter- 
ing (W). The crucial physical issues are revealed by the 
consideration of a ferromagnetic state in which all impu- 
rity spins Si are aligned, say, in the z direction. Then the 
carriers with spin parallel to Si feel a potential — J + W 
on each magnetic impurity site and anti-parallel carriers 
feel a potential J+W. These potentials self-consistently 
rearrange the electronic structure. The spin-dependent 
part of this rearrangement provides the energy gain which 
stabilizes the ferromagnetic state. The key physics issue 
is, evidently, whether the potential W + J is weak (so 
its effect on carriers near the lower band edge is simply 
a scattering phase shift) or strong (so only majority spin 
or perhaps both species of carriers are confined into spin- 
polarized impurity bands). Recent density functional su- 
percell calculations^ suggest that in GaMnAs — J + W 
is close to the critical value for bound state formation 
for the majority spin systems. Unfortunately, the precise 
effective values of J and W are typically unknown in a 
DMS system, and may have to be extracted experimen- 
tally. 

We assume that magnetic impurities under considera- 
tion enter substitutionally at the cation sites (e. g. Mn 
impurities at Ga sites) and the IIIi-^Mn^V system as 
a lattice of sites, which are randomly nonmagnetic (with 
probability 1— x) or magnetic (with probability x), where 
x now indicates the relative concentration (i.e. per Ga 
site) of active Mn local moments in IIIMnV. If more com- 
plete information about Mn locations on the GaAs lattice 
becomes available it will be straightforward to incorpo- 
rate that in the DMFT formalism. 

We now introduce the DMFT for the Hamiltonian 
given by Eq. (JTJ. Within the general scheme of the 
DMFT, the local (momentum independent) self energy 



of the system, E(ia>„), can be obtained from the time 
dependent mean field function. Then the single particle 
Green function is approximated by 



G(k, iu) v 



iuj n - (e k - fi) - E(iw n ) : 



(5) 



where fj, is the chemical potential. With the local self 
energy all of the relevant physics may be determined from 
the local (momentum-integrated) Green function defined 
by 



Gioc(«w n ) 



d 3 k 



G(k,iu n ) 



[ deD(e) 

J iuj n + [i- e 



- e - E a (iuj n ) 



(6) 



where D(e) = J d d p/ (2ir) d 6(e — e p ) is the density of states 
(DOS) for the noninteracting system. The information 
of the lattice geometry is included through the noninter- 
acting DOS. 

In our model G\ oc is a matrix in spin index and depends 
on whether one is considering a magnetic (a) or non- 
magnetic (b) site. Since G\ oc is a local function, it is 
the solution of a local problem specified by a mean-field 
function g , which is related to the partition function 
Z\ oc = J dS exp(— S'loc) with action 

Sloe = fdrj dr'J2ci(T)[g a 0af3 (T-T')]c p (r') 



a0 



f rfr^c+(r) [ JS • tT aP + W6 a> 



cp{r) (7) 



on the a (magnetic) site and 



S, 



loc 



dr / dr 



al3 



9Lp(T-T')ci(T)c p (T>), (8) 



on the non-magnetic (b) site. Here c a (r) (cj(r)) is the 
destruction (creation) operator of a fermion in the spin 
state a and at time r. <?o( T — T ') plays the role of the 
Weiss mean field (bare Green function for the local effec- 
tive action S'loc) and is a function of time. Gi oc depends 
only on frequency and is therefore the solution of a single- 
site problem. The local Green function G; oc of the effec- 
tive single-site problem is solely determined by the par- 
tition function, Z\ oc , namely, Gi oc (iuj n ) = 5 hi Z\ oc /5gQ 
which is identical to the local Green function computed 
by performing the momentum integral using the same 
self energy. Then the self-energy is defined by 



G lol, a f}( iuJ n)- 



(9) 



The a-site mean-field function g^ can be written as 
Soa/3 = ao+aiih-o" Q( g with m the magnetization direction 
and a\ vanishing in the paramagnetic state (T > T c ). 
Since the spin axis is chosen parallel to m g§ becomes a 
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diagonal matrix with components parallel (g^ = ao + ai) 
and antiparallel (g^ = a — a\) to rh. 

The form of the dispersion given in full Hamiltonian 
Eq. (1) applies only near the band edges. It is nec- 
essary for the method to impose a momentum cutoff, 
arising physically from the carrier band- width. We im- 
pose the cutoff by assuming a semicircular density of 
states D(e) = a% f j^jS(s - e pa ) = %/4i 2 - e 2 /2nt with 

t = (2-7T) 2 / 3 /ttioq. The parameter t is chosen to correctly 
reproduce the band edge density of states. Other choices 
of upper cutoff would lead to numerically similar results. 
This choice of cutoff corresponds to a Bethe lattice in in- 
finite dimensions. Other (perhaps more realistic) choices 
for the density of states would give results qualitatively 
similar to our results since the band edge density of states 
has the correct physical behavior in our model. For this 
D(e) the self consistent equation for go obeys the equa- 
tion 



lUJ n 



(1 



xXg^)- 1 ) 
g(iw„) + (j§-«r a /s+W 



(10) 



where the angular brackets denote averages performed 
in the ensemble defined by the appropriate Z\ oc , i.e., 
(A) = J dSP(Q)A, with P(Q) = exp(-5 loc )/Zi oc . With 
these mean field functions the local Green function can 
be written as 



G\ oc {iui n ) = iu n + (i - g . 



(11) 



The self energies are evaluated using Eq. and the 
full Green function from Eq. (JSJ). Physical observables 
can be obtained from the full Green function G(k, lo). 
In particular, the mean field function g§ a can be easily 
calculated at T = 



fftkrM =Ld + (1- X 

and T>T r 



gov - (W =f J) goa 



g 0tT (uj) = uj + ^-x- —2 (l-x) . 13 

(go„ - Wy - J 1 g a 



III. DENSITY OF STATES 

The density of states plays crucial roles in determining 
the physical properties of the DMS system. Especially, 
the formation of the impurity band arising from the im- 
purity doping gives rise to many different aspects from 
the continuum (i.e. virtual crystal approximation) semi- 
conductor band model. In this section, we calculate the 
DOS for different parameters (J, x, W, and T) and show 
how the impurity bands are formed and separated from 
the main band. We describe the DOS of dynamical mean- 
field calculations applying to simple semi-circle models. 
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FIG. 1: Temperature dependence of the DOS for x — 0.05 
and for a fixed coupling (a) J = l.Ot and (b) J = 2.0t. The 
evolutions of majority (minority) spin DOS show in top (bot- 
tom) panels for various temperature T/T c —0.0, 0.5, 0.7, 0.9, 
and 1.0. 



The DOS is given by the imaginary part of the Green 
function 



7T 



(14) 



In our model for a strong magnetic coupling J > J c 
two spin polarized impurity bands appear at the bottom 
(majority spin) and at the top (minority spin) of the 
main band. Each isolated impurity band has the weight 
x. However, if the coupling is not strong (J < J c ) the 
impurity band is not completely separated from the main 
bend. All DMS samples show that the carrier density 
is much smaller than the impurity concentration (n < 
x) due to the heavy compensation. As the system is 
the partially compensated, the chemical potential \l is 
located in the lower impurity band (if the impurity bands 
are formed) or the lower band edge (if the bands are not 
formed). Thus all physical properties arc determined in 
the lower energy band edge. Throughout this paper we 
only show the DOS near the lower energy band. 

In Fig. ^ we show the calculated DOS for various 
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FIG. 2: The calculated majority spin DOS at T = (upper 
half) and T = T c (lower half) is shown for various doping x = 
0.001, 0.01, 0.05, 0.1 and a fixed coupling (a) J = l.Ot, and 
(b) J = 1.51 



temperatures as a function of energy. The evolutions 
of majority (minority) spin DOS are shown in top (bot- 
tom) panels. In Fig. ^a) the strength of the coupling 
constant (J = l.Ot) is not strong enough to form the 
impurity band. Note that this value of coupling con- 
stant (J = l.Ot) is the critical value for impurity band 
formation as x — > 0. At T — the majority (minor- 
ity) spin band is shifted to lower (higher) energy com- 
pared with the noninteracting band which has a band 
edge at u> = —2.01 Thus all carriers are fully polarized 
when the carrier density is low. As the carrier density 
increases they start to occupy the minority band if the 
chemical potential crosses into the minority spin band. 
Recently we showed that the optical conductivity of the 
system is dramatically changed with the occupation of 
the minority bandit. We show in this paper that the 
calculated dc transport properties are also very sensitive 
to the minority band occupation. As the temperature 
increases, the minority band occupation grows and the 
carriers with minority spin increase due to thermal fluc- 
tuations. At T = T c both spins are equally populated 
and the bands becomes symmetric. As expected we have 
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FIG. 3: The calculated majority spin DOS at (a) T = 
and (b) T — T c is shown for a fixed x — 0.05 and for various 
coupling constant J/t — 0.5, 1.0, 1.25, 1.5, 2.0. 



a separated impurity band from main band for a strong 
coupling J = 2.0t shown in Fig. Q^b). When n = x the 
impurity band is fully occupied, no low energy hopping 
processes to main band are allowed and the system be- 
comes a band insulator. If the impurity band is partially 
occupied the derealization energy increases. At the half 
filling of the impurity band the system has the highest 
T^k. As the temperature increases spin disorder grows 
and the band becomes symmetric. But the impurity band 
width of the paramagnetic state (T > T c ) is smaller than 
that of the ferromagnetic state. This band shrinking oc- 
curs because the neighboring spins in the paramagnetic 
state are uncorrelated. 

In Fig. |3we show the calculated majority spin DOS 
corresponding to the disordered spin state (at T = T c , 
bottom panels) and ferromagnetic state (at T = 0, top 
panels). The evolutions of the energy (uj) dependent DOS 
are shown for different doping parameter x =0.001, 0.01, 
0.05, 0.1 and for fixed coupling constant (a) J = l.Ot and 
(b) J = 1.51 In our model the impurity level (acceptor 
energy level) and the formation of an impurity band de- 
pend on the ferromagnetic coupling J. If J < J c = t the 
impurity level is not isolated from the main band, but if 
J > J c we find an isolated impurity level below the main 
band. The small dot in Fig. |2Ib) indicates the isolated 
impurity level in the dilute limit (x — > 0). As a; increases 
for J > J c an impurity band centered around the impu- 
rity level is formed below the main band. For J < J c 
the impurities give rise to band tailing in the main band 
edge instead of forming impurity band. The band width 
of the impurity band for J > J c increases with x since 
the number of states of the impurity band increases with 
x. If x is bigger than x c the impurity band merges into 
the main band. For J — l.ht we have x c = 0.032 at 
T = OK and x c = 0.071 at T = T c . 
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0.00 0.02 0.04 0.06 0.08 
FIG. 4: The majority (top panel) and minority (bottom H 
panel) DOS at T = are shown for a fixed value of J = 1.5t 

and x = 0.05 for various potential scattering W/t = -0.5, FIG - 5: The density dependent resistivity at T = T c for 
—0.25 0.0 0.25 0.5. a fixed value of x — 0.05 and for various coupling constant 

J/t = 0.7, 1.0, 1.2, 1.4, 1.5, 2.0. Insets show the DOS with 
the corresponding parameters. 



In Fig. [3]we show the majority spin DOS at (a) T = 
and (b) T = T c for a fixed x = 0.05 and for various 
coupling constant J/t = 0.5, 1.0, 1.25, 1.5t, 2.0t. For 
J < J c we see the expected band shift proportional to xJ . 
For J > J c an impurity band centered at impurity level 
and containing x states is seen to split off from the main 
band. Due to the widening of the impurity band with x 
we find the separated impurity band when J > 1.58i at 
T = and J > 1.45* at T = T c for x = 0.05. Thus the 
separated impurity band is expected with small x and 
large values of J. 

In Fig. 0|we show the calculated DOS at T = for a 
fixed value of J = l.bt and x = 0.05, and for various val- 
ues of potential scattering W =—0.5t, — 0.25t, 0.0, 0.25, 
0.5. The DOS of majority (minority) spin is shown in 
the top (bottom) panel. When we include the potential 
scattering (W) in addition to the spin-spin coupling (J), 
the carriers with spin parallel to local impurity spin feel 
a potential —J + W on each magnetic impurity site and 
anti-parallel carriers feel a potential J + W . Thus the 
formation of the impurity band and the corresponding 
physical properties of the system depend on the com- 
bined coupling W ± J. Fig. 0] shows that while the band 
edge of the minority spins is slightly dependent on the 
potential scattering, the majority spin DOS is strongly 
affected by the potential scattering. Even weak potential 
scattering can change the extended majority spin band 
into the spin-polarized impurity bands or the well formed 
impurity band into the band tail of the main band. 

In the following section we show that our calculated 
resistivities depend strongly on whether the carriers are 
within the impurity band or in the band tail of the ex- 
tended main band. 



IV. DC-RESISTIVITY 

The conductivity is calculated from the usual Kubo 
formula. The Kubo formula for the conductivity a in- 
volves the two particle current-current response function. 
Since the irreducible vertex in the response function is 
purely local in our approximation of DMFT there is no 
vertex correction 18 . Thus, within this approximation 
only the simple bubble survives and the real part of the 
finite frequency conductivity is given by 

MT) - / AJM .(.)/±i ^»-g" + n " 

x A tT (e,uj)A <T (e,uj + fl), (15) 

where A a (e 7 uj) = — (l/7r)ImG cr (e, lo) is the spectral func- 
tion, (f>(e) = (At 2 — £ 2 )/3 is the current vertex for the 
Bethe latticeiS,, f(u>) is the Fermi distribution function. 
For the hypercubic lattice^ 3>(e) = 1 has been used in 
Eq. (JTSJ, but for the Bethe lattice the explicit form 
is derived in ref. yjj The dc-resistivity p = l/a^c 
is then found from Eq. Ijl5|l in the limit O — ^ with 
Vdc = cr(Q — > 0), which is given by 

^c(T) = e 2 £| deD(eMs) J ^ ^( e>W ). 

(16) 

In Fig. 03 we show the calculated resistivity at T = T c 
as a function of carrier density for a fixed value of 
x = 0.05 and various coupling constants J/<=0.7, 1.0, 
1.2, 1.4, 1.5 and 2.0. Insets show the density of states 
near the band edge corresponding to the disordered spin 
state (T — T c ). All calculated results show that p di- 
verges as n -> due to the absence of carriers. As den- 
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FIG. 6: The calculated resistivity for (a) a strong coupling 
J = 2.0t with densities n = 0.01, 0.015, 0.02, 0.025, 0.03 
0.035, 0.04, 0.045 (top to bottom)and (b) a weak coupling 
J = 2.0t with densities n = 0.008, 0.01, 0.015, 0.02, 0.025, 
0.035, 0.04, 0.045 (top to bottom). Dots indicate T c for given 
densities. 



sity increases we find two different behaviors depending 
on the formation of the impurity band. (The critical 
coupling constant which gives rise to the formation of 
the well separated impurity band below the main band 
for x = 0.05, T = T c is J c = 1.48t.) In the weak cou- 
pling limit J < l.Ot, where the impurity band formation 
dose not happen, the resistivity decreases as the density 
increases monotonically. However, in the strong coupling 
limit (J > J c ), where the impurity band is formed, the 
resistivity diverges again when the impurity bands are 
fully filled (i.e., for n = x) and the system becomes 
a band insulator. In the intermediate coupling regime 
(l.Ot < J < J c ) the resistivity shows non-monotonic be- 
havior. 

In Fig. we show the calculated resistivity as a func- 
tion of temperature for various density. For the strong 
coupling limit (J = 2.0t), in Fig. Eta), we find a crossover 
(resistivity peak) separating a good metal at low T from a 
semiconductor at higher T. The resistivity peaks are pro- 



portional to the energy separation between the chemical 
potential and the band edge of the minority band. In the 
high temperature regime above the resistivity peak the 
decreasing resistance with increasing temperature, char- 
acteristic of a semiconductor or insulator, is due to ther- 
mal excitation of the carriers from impurity band to the 
upper minority spin band. The resistivity decreases with 
density because the carriers in the main band are scat- 
tering off the impurities. But at low temperatures most 
carriers in the impurity band contribute to the scattering 
and give non-monotonic density dependence of the resis- 
tivity. The metallic behavior at low temperature can be 
understood by the disappearance of the coherent central 
quasiparticle peak in the DOS. For the weak coupling 
limit (J = l.Ot), in Fig. HJb), the resistivity can be ex- 
plained by the scattering effects in the main band except 
the behavior in the ferromagnetic state (T < T c ). For 
T <T C the resistivity changes from insulating at low den- 
sities to metallic at high density. Details of this behavior 
are given in Fig. [7| We see that the resistivity crossover 
takes place only at low densities since the minority band 
is occupied by the carriers at high density. 

In Fig. [7| we show the dc-resistivity as a function of 
temperature for different densities in the low tempera- 
ture regime (T < 2T C ). In Fig. Eta) we use the param- 
eters x — 0.05 and a strong coupling J = 2.0t. In this 
case all carriers (if n < x) occupy the impurity band 
and stay mainly at Mn sites. Thus, the carriers in the 
impurity band follow the fluctuation of the localized Mn 
spin. At half filling of the spin polarized impurity band 
(n = 0.025) the resistivity has the lowest value and its be- 
havior corresponds to that of the double exchange (DE) 
model (see Fig. ^J|, that is, the resistivity decreases be- 
low the critical temperature because of the spin disorder 
scattering as the temperature decreases^. In the low 
density limit (n < 0.25, less than half filling) the resis- 
tivity is dominated by the 'impurity band' contribution 
and the resistivity increases as density decreases due to 
the lack of mobile carriers. As the carrier density is in- 
creased above the half filling (n > 0.25) the resistivity 
increases due to the filling of the band. We also find 
a very different temperature behavior of the resistivity 
from the low density case. As the temperature decreases 
the resistivity increases just below T c , then decreases at 
very low temperatures. The counter- intuitive increase of 
resistivity just below T c (as T is decreased) arises be- 
cause, as the carrier spins are aligned to the impurity 
spins, the binding of the carriers to the impurity spins 
increases corresponding to an increase in the basic scat- 
tering rate. In the very low temperature range, however, 
the DE-like mechanism dominates, which gives rise to 
decrease of the resistivity. (When the impurity band is 
spin-polarized, carriers which are bound to impurity site 
must have spins parallel to impurity spin. Thus, as the 
spins order ferromagnetically, the basic ability of carri- 
ers to move in the impurity band increases.) The overall 
temperature dependence of resistivity is very weak in the 
strong coupling limit. This weak T-dependence of the dc 
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FIG. 7: The low temperature dc-resistivity as a function of 
temperature for a fixed x = 0.05 and for various densities. 
In (a) the results for a strong coupling constant (J = 2.0f) 
and densities n = 0.01, 0.012, 0.015, 0.025, 0.3, 0.35, 0.38, 
0.40 are shown. In (b) we use J = l.Ot and n = 0.015, 0.016, 
0.0163, 0.0165, 0.18, 0.2, 0.25, 0.3 (from top to bottom) . Dots 
indicate the T c for given parameters. 



resistivity below T c occurs because the increase in scat- 
tering rate due to the binding is compensated by DE-likc 
mechanism. 

In Fig. |7Jb) we show the results for a weak coupling 
limit, J = l.Ot. In this case the impurities contribute 
to form the band tail of the main band and the chemi- 
cal potential lies in the main band. In the low density 
limit (n < n c = 0.0164, where n c is the density above 
which the minority band starts to fill) the resistivity be- 
low T c increases as the temperature decreases. The in- 
crease in resistivity is due to increased carrier-spin cou- 
pling as mentioned above. But, in the high density limit 
(n > n c , where at T — the minority spin-band is oc- 
cupied) we find the resistivity decreasing as the temper- 
ature decreases. This metallic behavior in high density 
regime and for a low coupling constant can be understood 
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FIG. 8: The energy evolution of optical conductivity for (a) 
a intermediate coupling J = l.Ot and (b) a strong coupling 
J = 2.0t. Solid (dashed) lines correspond to the results for 
T = (T — T c ). Insets show the dc conductivities with 
corresponding parameters. Arrows indicate T c . 



by the small scattering rate of the carriers in the minor- 
ity band. If the carriers are near the edge of the majority 
spin band, the carriers form a spin-polarized bound state, 
so the effective scattering rate strongly increases, which 
gives rise to the insulating behavior as shown in the low 
density results. On the other hand, if the carriers are in 
the minority spin band the carriers form an anti-bound 
state at the top of the band, so that at the physically 
relevant lower band edge, the effective scattering rate de- 
creases. In addition in three dimensions the vanishing of 
the density of states at the band edge further decreases 
the scattering. These effects are quite large and domi- 
nate as the minority band is occupied, which gives rise 
to the metallic behavior in the high density limits. In 
the paramagnetic state the resistivity is almost temper- 
ature independent. In the limit T > T c the local Green 
function G(u) is temperature independent and the resis- 
tivity depends on temperature only weakly through the 
Fermi function (thermal smearing around the chemical 
potential) . 

In Fig. [SJwe show the relation between the optical con- 
ductivity and dc-resistivity. The main panels of Fig. |S| 
show the evolution of the conductivity for two couplings; 
weak (J = l.Ot, Fig. Ufa), where the impurity band for- 
mation is not accomplished), and strong (J = 2.0t, Fig. 
|SJb), where the impurity band is well formed); the insets 
show the dc resistivity with the same parameters. The 



9 



j=t x=0.05 n=0.015 



28.6 



28.2 



17.0 
16.0 

15.0 

14.0 
9.0 



H 



8.0 - 



7.0 
6.5 



5.5 

4.5 
3:0 

2.0 - 



1.0 
0.00 





W=-0.5t ^ 

t 




W=-0.2t 

• 


• 

• 

•« 


W=0.0 

J : 




J «. 

/ t 




w=o.o5t ; 




/! 


W=0.2t 



0.005 



0.01 

T/t 



0.015 



0.02 
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FIG. 10: Resistivity p(T, B) for various magnetic field is 
given as a function of temperature for (a) J = 2.5i, x = 1.0, 
n = 0.5, and (b) J = 2.5i, x = 0.05, n = 0.025. Insets 
show the normalized magnetization M/Mq, where Mo is the 
saturation magnetization. 



solid (dashed) lines indicate the optical conductivity for 
T = (T = T c ). When the impurity band is not formed 
(for J = l.Ot) we find approximately the Drude form for 
optical conductivity expected for carriers scattering off 
random impurities (a closer examination reveals minor 
differences due to density of states variations near the 
band edge). In this case (and for n < n c ) the dc conduc- 
tivity shows insulating behavior due to the formation of 
the bound state. Since the carrier density is low enough 
not to fill the minority band, the formation of anti-bound 
state dose not take place, which reduce the scattering 
rate. In the J — 2.0t case the density of states plot shows 
the formation of an impurity band and the corresponding 
conductivity has two structures: a low- frequency quasi- 
Drude peak corresponding to motion within the impurity 
band and a higher frequency peak corresponding to exci- 
tations from the impurity band to the main band. In this 
case the dc resistivity shows metallic behavior because 
the reduction of spin-disorder scattering dominates over 
the bound-state formation. In Ini^Mn^As^ we find the 
Drude-like conductivity is correlated with the insulating 
behavior, but in Gai-ajMn^ the mid infrared peak in 
the optical conductivity is closely related to the decrease 
of the dc-resistivity below the critical temperature. 



Fig. shows the sensitivity of the predicted behavior 
to potential scattering. In this figure we use the parame- 
ters: x — 0.05, n — 0.01, J — l.Ot, and for various poten- 
tial scatterings W/t =-0.5, -0.2, 0.0, 0.05, 0.2. At zero 
scalar potential (W — 0.0, middle panel) the impurity 
band is not formed, and the resistivity shows insulating 
behavior below T c due to the bound state formation of 
the carrier spins with impurity spins. As the potential 
is made more attractive (negative), the impurity band 
features become pronounced and the spin-disorder scat- 
tering decreases (the carriers move easily in the impurity 
band). As the potential is made more repulsive (posi- 
tive), the impurity band rapidly rejoins the main band. 
This reduces the energy gap between the Fermi energy 
and the band edge of the minority band, and for large 
enough repulsive potential the minority band starts to be 
filled by the carriers. Thus, the anti-bound state forma- 
tion dominates and it reduces the scattering rate below 
T c and shows a metallic behavior. The DMS transport 
properties are sensitive to the combined coupling W ± J, 
and not solely to the exchange coupling J. 

Now we compare the transport properties of DMS to 
those of another system with strong carrier-spin cou- 
plings, namely CMR manganites 10 . In CMR, instead 
of being dilute random impurities as in DMS, the Mn 
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ions form an ordered lattice. They possess large local 
moment, to which mobile carriers are very strongly cou- 
pled. Thus instead of a spin-polarized impurity band, 
there is a spin-polarized conduction band, sufficiently 
well separated from other spin bands. The periodic ar- 
rangement of the Mn sites means that (in the absence 
of other physics) the scattering rate decreases as T is 
lowered. CMR materials can be understood well by the 
double exchange (DE) models. In our model this cor- 
responds to x = 1, in which all magnetic ions replace 
Ga at the cation sites. In this case (x — 1) the fully 
polarized spin band is well separated from the other 
bands instead of forming impurity band. The temper- 
ature dependences of the resistivity p(T, B) and magne- 
tization M(T, B) are given in Fig. HOt a) for x = 1.0 with 
J = 2.5t, n = 0.5, and various magnetic fields. Above T c 
the resistivity has a small temperature dependence since 
the local spin fluctuation is saturated above T c . Below 
T c resistivity decrease as magnetization increases. The 
origin of the resistivity dependence on the magnetization 
is spin disorder scattering, which gives rise to the scaled 
behavior of the resistivity, p(M)/p(M = 0) = 1 - CM 2 , 
where C is a temperature/field independent constant 20 . 
The origin of the resistivity is qualitatively explained by 
the carrier scattering due to the thermally fluctuating 
spin configurations, or the spin-disorder scattering. As 
the spontaneous or the induced magnetic moment is de- 
veloped, the amplitude of the spin fluctuation decreases 
so that the resistivity also decreases. In Fig. HOf b) we 
show the resistivity for DMS system with x — 0.05. The 
overall features of temperature and magnetic field depen- 
dence look similar to the DE model. However, we find 
that the negative magnetoresistence at T c is very weak 
and the resistivity above T c is not saturated. The fast 
drop of the resistivity just below T c can be explained by 
spin-disorder scattering. These ideas of DE model have 
limited applicability to the DMS systems, i.e. only for 
strong couplings and near half filling of the spin polar- 
ized impurity band. As shown in previous figures these 
ideas cannot explain the resistivity behavior in low cou- 
pling and high density regimes. Note that an important 
ingredient of DMS transport properties is missing from 
our DMFT theory which, while accounting well for the 
non-perturbative effects of spin disorder (</) and local po- 
tential (W) scattering by the magnetic impurities, leaves 
out all ionized impurity disorder that may very well be 
important. 



V. ELECTRICAL RESISTIVITY: BOLTZMANN 
TRANSPORT APPROACH 

When the dominant scattering mechanism is the scat- 
tering by charged impurities the Boltzmann transport 
theory may be used to calculate the electrical resistiv- 
ity of the carriers since DMFT is not well-suited for 
treating long-range disorder. Due to the band splitting 
in the ferromagnetic state the carrier densities n± for 



spin up/down are not equal. Note that the total density 
n = n + +rt_. In this situation the total conductivity can 
be expressed as a sum of spin up/down contributions 

a = (T + + er_, (17) 

where a± is the conductivity of the (±) spin subband. 
The conductivities a± are given by 



n±e'{r±) 



(18) 



where m is the carrier effective mass, and the energy aver- 
aged transport relaxation time (r±) for the (±) subbands 
are given in the Boltzmann theory by 



(r±) 



J deT±(s)e 


de 


J dee 


\ <9/±( £ )~ 
de 





(19) 



where T±(s) is the energy dependent relaxation time for 
the (±) subbands, and / ± (e) is the carrier (Fermi) dis- 
tribution function 



1 



1 + e [£-p±(T)]/fc B T> 



(20) 



where fx (T) is the chemical potential at finite temper- 
ature. The energy dependent relaxation time is given in 
the Born approximation by 

MeiOP 1 = Y J ^mV^\ 2 (l~cos9 kk/ )6(e k -e kl ), 

(21) 

where A,; is the charged impurity concentration, and 
Vk-k' i s the screened carrier-impurity Coulomb interac- 
tion, which can be expressed as 



V sc (q) 



4nZe 2 



1 



Kq 2 l + [q s {q)/q] 



(22) 



where Z is the charge of impurities, n the background 
(GaAs) lattice dielectric constant, and q s the tempera- 
ture dependent screening function. 

In this paper we consider two screening approaches: 
Thomas-Fermi (TF) approximation and random phase 
approximation (RPA). Within TF screening we have 



Itf 



k° F 



k° F 



(23) 



where q^ F — ikp/iras with k F = (^TT 2 n) 1 ^ 3 being the 
Fermi wave vector of the spin unpolarized state and cle = 
nh 2 /me 2 the effective Bohr radius. k F = { t oTT 2 n±) 1 ^ is 
the Fermi wave vector of the each spin-split subband. 
Note that TF screening wave vector, being a long wave- 
length approximation, is dependent only on the spin 
polarization, but not on the temperature explicitly. It 
is therefore temperature independent above the critical 
temperature. Within TF screening approximation we 
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FIG. 11: Wave vector dependence of the effective screening 
function (polarizability function 11(g) normalized its param- 
agnetic state value N F = II para (0) = mk^/n 2 ) for various 
polarizations 5 — (n+ — n_)/n. 



have the energy dependent relaxation time (by integrat- 
ing Eq. 21) 



3 n h 



k F 



(24) 



where k F = (k+ + k F )/2 and F(x) = (2/s 2 )[ln(l + x) 
x/ + (l 
kp = 



x)]. When a system is fully polarized we have 
l / 3 k F and q s = q F F /2^^ 3 . For the random phase 
approximation (RPA) the temperature dependent screen 
function is given by 



Qtf 



^U + ( q> T) + ^U.( q ,T) 

hjp h,p 



, (25) 



where TL±(q,T) is the temperature dependent static 
Lindhard function for each spin subband. In Fig. 1111 
we show the screening function q 2 (q,T — 0)/q TF = 
U(q)/N F , where N F = IP ara (0) = rnkp/ir 2 , for vari- 
ous polarization 5 = {n + — nJ)/n. For an unpolarizcd 
state (5 = 0) there is an inflection point at q = 2k F , 
which is the usual Kohn anomaly. For partially polar- 
ized state < 5 < 1 we have two inflection points at 
q = 2kf = (2n ± /n) 1 / 3 and the value of U(q = 0)/N F 
decreases as S increases. When the system is fully polar- 
ized 6 = 1, and 11(g) has an inflection point at q = 2k F , 
and U(q = 0) = N F /2 2 / 3 . 

At T = 0, f ± (e) = 9(Ep - e) where E$ is the Fermi 
energy for the (±) subbands, and then (t±) ee t(E f ) 
giving the familiar result a ee p^ 1 = n + e 2 T(E F )/m + 
n^e 2 T(E F )/m. Within TF screening the ratio of the 
resistivity of the fully polarized state to that of the 
unpolarized state becomes p(5 = l)/p(S = 0) = 
2 5 / 3 F(2 4 / 3 a; )/F(a; ), where x = {2k F /q TF ) 2 . (Note 
xq oc n 1 / 3, .) In Fig. ^| we show the calculated re- 
sistivity as a function of 5 = {n + — nJ)/n for sev- 
eral value of xq. Solid (dashed) lines show the results 
calculated with RPA (TF) screening. In general, for 
small values of xq, p(5 = 1) > p(5 = 0), but for 




FIG. 12: Calculated resistivity as a function of S = (n + — 
n-)/n for various value of xo = 1.0, 2.0, 5.0, 10.0 (top to 
bottom). Solid (dashed) lines show the results calculated with 
RPA (TF) screening function. 



1) < p(S = 0). At xq w 3 we have 
= 0). For GaAs xq = 3 corresponds to 



large x , p{8 = 
p(S = l)=p(S 
the hole density n w 10 19 cm 3 . In relatively high den- 
sity limits (large xq) the two approximations agree very 
well, which indicates that the q = scattering mostly 
contributes to the scattering time. However, in the low 
density regime (small xq) we that find the two screening 
theories give very different results for the spin polarized 
state. Noting that the TF approximation is just the long 
wavelength limit of RPA, we emphasize that RPA is ob- 
viously the better theory. 

In Figs. El an d E] we show our calculated resis- 
tivity for GaMnAs samples. We use the parameters 
corresponding to GaMnAs: dielectric constant k = 
12.9, hole effective mass m = 0.5m e , impurity density 
N{ = 10 21 cm -3 , and the magnetic coupling constant 
J =120 meV nm 3 . Fig. 1131 shows normalized resistivity 
p(T) I p(T c ) as a function of temperature T/T c for several 
values of hole density. In Fig. I13l there are two indepen- 
dent sources of temperature dependence in our calculated 
resistivity — one source is the energy averaging defined in 
Eq. (|19fl and the other is the explicit temperature depen- 
dence of the finite temperature RPA screening function 
q s (q, T) . Since the Fermi temperature is much higher 
than the critical temperature (T C /T F -C 1) the screen- 
ing function is weakly temperature dependent above the 
critical temperature (unpolarized state). Thus the de- 
crease of the resistivity above T c with increasing temper- 
ature arises from the thermal energy averaging. This is a 
well-known high-T effect of ionized impurity scattering in 
semiconductors: p(T) decreases with increasing T due to 
the thermal averaging. As density decreases this effect is 
enhanced. In the ferromagnetic polarized state (T < T c ), 
as the temperature decreases from T = T c , the screen- 
ing length increases until all the holes are polarized. In 
this temperature range we find a strong temperature de- 
pendent resistivity. This strong temperature dependence 
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FIG. 13: Calculated resistivity, p(T)/p(T c ), as a function 
of temperature. In (a) the RPA results are shown with hole 
densities p/Ni = 0.05, 0.1, 0.2, 0.5. In (b) the TF screening 
results are shown with hole densities p/N = 0.01, 0.05, 0.1, 
0.2. Here the impurity density Ni — 10 21 cm' 3 is used. 



arises from the low temperature screening function and 
the change of the carrier densities of the each subband. 
When all the carriers are polarized (i.e. at low density 
and temperature) the screening function is almost inde- 
pendent of temperature and the resistivity is also there- 
fore temperature independent. The interplay between 
the screening length and the down spin (-) carrier den- 
sity gives rise to the minimum of the resistivity in the low 
density regime. The decrease of the scattering time due 
to the fast increase of the down spin density overwhelms 
the increase of the scattering time due to the decrease of 
the screening length, producing the non-monotonic be- 
havior for p/N t < 0.1 and T/T c ~ 0.5 in Fig. 13. But 
in the high density regime, where the spins are not fully 
polarized even at low T, screening is the dominant ef- 
fect on the temperature dependent resistivity. As the 
density increases the relative low temperature resistivity, 
p(T)/p(T c ), decreases until the holes are partially polar- 
ized. The change of screening wave vector is larger in 
this case, leading to the monotonically increasing p(T) 
in the T < T c regime. 

In Fig. ^]we show our calculated resistivity as a func- 
tion of temperature for two hole densities (a) p/N = 0.1 
and (b) p/Ni —0.2 with finite temperature RPA screening 
function. The temperature dependence of the impurity 
(i.e. Mn) moment magnetization as well as the hole spin 
polarization is given in the insets, which is calculated us- 
ing the Weiss molecular mean-field theory for delocalized 
carriers^. Note that the hole gas is almost fully spin po- 
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FIG. 14: Calculated resistivity as a function of temperature 
for a hole density (a) p/rii = 0.1 and (b) p/Ni —0.2 with 
finite temperature RPA screening function. The top insets 
show the temperature dependence of the Mn ions and hole 
magnetization. The bottom insets show the calculated tem- 
perature dependence of the TF screening wave vector g s (0). 
Here the impurity density Ni = 10 21 cm -3 is used. 



larized upto T/T c = 0.5 at low density, but at higher 
density the holes are partially polarized even at T = 0. 
The bottom insets show the calculated temperature de- 
pendence of the finite temperature screening wave vector 
g s (0). At high density the screening wave vector changes 
by 6% when the temperature increases from zero to T c 
due to the partial polarization of the holes at T = 0. But 
at low density the decrease of the screening wave vector is 
about 20%. In the metallic GaMnAs samples the change 
of the resistivity when the temperature goes from T c to 
zero is about 20% in good agreement with our calcula- 
tion. Similarly the observed decrease of p(T) for T > T c 
also arises naturally in our theory as a consequence of 
thermal averaging. Thus, the temperature dependence 
of the resistivity in the metallic DMS samples may be 
arising almost entirely from the temperature dependence 
of screening and thermal averaging in the charged impu- 
rity scattering. 



VI. DISCUSSION 

We have developed in this work two complementary 
theories for understanding DMS transport properties. 
Our work establishes that DMS transport, even in ideal 
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intrinsic circumstances, is rather complex, and depends 
sensitively on many system parameters including the ex- 
change coupling, the magnetic impurity density, the car- 
rier density, the temperature, the band parameters of the 
parent semiconductor (e.g. effective mass, etc.), and the 
details of the charged impurity distribution (and there- 
fore compensation). Given such a complex parameter 
space, it is not meaningful to try to develop quantita- 
tive transport theories at this early stage of the sub- 
ject since all the intrinsic parameters may not be known. 
We therefore focus in this work on developing a compre- 
hensive qualitative theoretical description which empha- 
sizes general broad features of how various parameters 
affect DMS transport behavior. As such we have concen- 
trated in this work on understanding temperature and 
carrier density dependence of dc resistivity in a model 
DMS system, keeping primarily the extensively studied 
Gai-^Mn^As system in mind. Even for GaMnAs, the 
transport data for various values of Mn concentration 
(x « 0.01 — 0.1) cover much too broad a range of behav- 
ior for a unified and coherent theoretical description. For 
example, low (and sometimes large) Mn concentrations 
(x < 0.03 and sometimes x > 0.05) are known to lead 
to insulating transport behavior usually attributed to a 
disorder-driven metal-insulator (Anderson) localization. 
We neglect all localization effects in our theory. The lo- 
calized GaMnAs regime in all likelihood requires its own 
characteristic theoretical description such as the bound 
magnetic polaron percolation theory^ whose transport 
properties^ have recently been theoretically analyzed. 

Even without the disorder-induced strong localization 
complications, neglected completely in this work, we face 
the formidable difficulty of using the semiconductor va- 
lence/conduction band (the valence band for GaMnAs, 
where the carriers mediating the ferromagnetic interac- 
tion are holes) or the magnetic impurity induced impu- 
rity band (i.e. Mn induced d-band in the fundamental 
band gap of GaAs) picture for describing the carrier dy- 
namics. The precise nature (i.e. valence band versus 
impurity band) of the DMS carriers is still a controver- 
sial issue although it is likely that at the very high doping 
densities (e.g. Mn density ~ 10 21 cm -3 in GaMnAs) of 
DMS interest the impurity band overlaps strongly with 
the valence band (i.e. forms the tail of the valence band), 
and therefore, the distinction between the valence and 
the impurity band picture is not a real qualitative dif- 
ference. Our DMFT theory, presented in sections II-IV 
of this paper, clearly shows that in the strong exchange 
coupling (J/t S> 1) regime the impurity band picture ap- 
plies whereas in the weak-coupling regime (J/t < 1), the 
valence band picture applies. It is possible that GaMnAs 
belongs to the intermediate coupling regime (J/t ~ 1), 
where it may be more appropriate to think of the holes 
to be residing in the extended tail of the valence band, 
presumably with an enhanced effective mass compared 
with the GaAs valence band hole mass. Such a coupled 
impurity-valence band picture of GaMnAs is consistent 
with recent optical spectroscopy measurements^!, but 



more experimental work is needed to settle this question. 

The theoretical strength of our DMFT description is 
that, being a nonperturbative technique, it can handle 
both strong-coupling and weak-coupling regimes, and our 
results presented in Figs. 6 - 10 of this paper show 
qualitative difference between the strong-coupling regime 
(J/t = 2) with an impurity band well-separated from 
the semiconductor band and the intermediate-coupling 
regime (J/t = 1) with only band-tailing and no sepa- 
rate impurity band formation. Temperature, carrier den- 
sity, and impurity concentration all play qualitatively im- 
portant roles in determining the dc transport properties 
within DMFT, and sorting out the details with respect 
to experimental results may be extremely difficult. 

The weakness of DMFT is that it can only include 
spin-disorder scattering (controlled by J) and Mn impu- 
rity induced local potential scattering (controlled by W) 
effects on transport properties. As such it leaves out the 
most important scattering mechanism which may be op- 
erational in real samples, namely scattering by charged 
impurities which is often the most important resistive 
scattering process in heavily doped semiconductos below 
the room temperature (or the optical phonon scattering 
regime which may well be above the room temperature) . 
The reason DMFT is unable to account for charged impu- 
rity scattering is that the Coulombic impurity potential 
is long ranged, and DMFT by construction is a local the- 
ory. Thus, rewriting our starting Hamiltonian (Eq. 1) 
more completely we have 

H = H host + H M + H AF + [Hi + H c ] , (26) 

where Hi is the carrier-charged impurity interaction and 
H c is the carrier-carrier (i.e. hole-hole in GaMnAs) in- 
teraction. In principle, the terms (i.e. H^, H c ) within 
the square bracket are parts of the Hhost, but it is im- 
portant to appreciate their considerable (perhaps even 
dominant) importance in determining the dc transport 
properties. To include the charged impurity scattering 
effects on transport, we use the highly successful and ro- 
bust semiclassical Boltzmann transport theory to DMS 
systems assuming a mean-field approach where the long- 
range Coulomb impurity potential arising from Hi (as- 
suming random impurity scattering) is screened by the 
polarization bubble diagrams arising from H c . This type 
of Boltzmann transport theory is extremely successful 
in describing semiconductor transport properties^. We 
note that our DMFT and Boltzmann transport theories 
are complementary - DMFT treats the spin disorder and 
the local potential scattering associated with Mn dopants 
and Boltzmann theory treats the scattering by screened 
charged impurity scattering. 

The magnetic DMS properties enter the Boltzmann 
theory only indirectly through the carrier spin polar- 
ization calculations. Spin disorder scattering is not ex- 
plicitly included in the Boltzmann theory although it is 
straightforward to do so. Our Boltzmann theory mani- 
fests nontrivial interplay among temperature dependent 
screening, temperature dependent spin polarization (i.e. 
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FIG. 15: Schematic diagram for the experimental GaMnAs 
dc-resistivity for decreasing hole density (from bottom to top) 
(see ref. [3]). 



spin up-down carrier densities), and thermal energy av- 
eraging, leading to temperature dependent resistivity 
(Figs. 13 and 14) that are rather similar to experimen- 
tal observations^ in GaMnAs. Based on this qualitative 
similarity we conclude that much of GaMnAs DMS trans- 
port is dominated by screened charged impurity scatter- 
ing with spin disorder scattering playing only a rather 
minor quantitative role. Recently, Lopez-Sancho and 
Brcy 15 have come to a very similar conclusion for GaM- 
nAs transport properties. 

Before concluding we now discuss our theoretical re- 
sults in light of the existing DMS transport data. Al- 
though there are some experimental transport data in 
other DMS systems (most notably InMnAs with qual- 
itatively similar behavior to GaMnAs), truly exten- 
sive and reliable transport data 3 are available only for 
Gai^Mn^As (in the x w 0.01 - 0.08 regime) DMS sam- 
ples. Even for this well-studied system, experimental 
transport results are problematic, and have considerable 
spread in the sense that nominally "identical" GaMnAs 
samples (i.e. same nominal carrier density and Mn con- 
centration made in the same growth run) may have dif- 
ferent T c and transport properties. This situation is im- 
proving as sample quality and processing (e.g. anneal- 
ing) techniques improve, but experimental DMS trans- 
port properties are still not robust in a quantitative way. 
This is of course very understandable given the very large 
parameter space (i.e. exchange coupling, hole density, 
Mn concentration, defects and impurities, compensation, 
band structure parameters, etc) that DMS transport de- 
pends on. With these serious caveats we show in Fig. 15 
a schematic depiction of the generic experimental obser- 
vation for p(T) in GaMnAs for various hole densities. At 
high hole density the system shows "metallic" behavior 
for T < T c with p(T) increasing somewhat with T upto 
T c , and then decreasing slowly for T > T c . Even the 
optimally doped most "metallic" GaMnAs is, however, 
at best a bad metal with mobilities of the order of 10 
cm 2 /Vs with kpl ~ 1 where I is the transport mean free 
path. It is important to realize that DMS transport is 
always extremely highly resistive due to the very large 



amount of impurities and defects invariably present in 
the low temperature MBE process needed for producing 
homogeneous GaMnAs samples. Thus, from the perspec- 
tive of heavily doped semiconductors, although the DMS 
systems may be above the Mott limit (i.e. the carrier 
density in the Mott metallic range), they are close to be- 
ing Anderson insulators due to strong disorder effects. As 
hole density decreases the system eventually becomes an 
insulator at low enough carrier densities (the top curve in 
Fig. 15), with p(T) decreasing monotonically increasing 
T. 

Two generic features of Fig. 15 are: (a) a peak (or a 
kink) in p(T) near T w T c , and (2) the slow decrease of 
p(T) for T > T c . Both of these generic features of p(T) as 
well as the "metallic" high-density behavior are qualita- 
tively (perhaps even semi-quantitatively) well-explained 
by our Boltzmann theory approach including only scat- 
tering by ionized impurity scattering. This is apparent by 
comparing Fig. 15 with Figs. 14 and 13 where our Boltz- 
mann theory results are shown. Physically, the increasing 
p(T) with T (< T c ) arises from the decreasing strength 
of screening due to the interplay of two independent and 
competing physical effects: Temperature induced sup- 
pression of screening and spin polarization (i.e. the spin 
polarization decreasing with increasing temperature) in- 
duced enhancement of screening with increasing temper- 
ature. As explained in Sec. V the competition between 
these two effects depend on the carrier density, leading to 
some weak non-monotonicity in p(T) for T < T c . In this 
screening picture, the resistivity peak or cusp at T w T c 
arises from the ferromagnetic to paramagnetic transition 
which affects screening as the carrier density of states 
(which is inversely proportional to the screening length) 
changes from one in the fully spin polarized state to two 
in the fully paramagnetic phase. Thus, even without any 
spin disorder scattering effects, just screened ionized im- 
purity scattering by itself will give rise to the peak or the 
kink in p(T) at T — T c . We believe that spin disorder 
scattering, which also produces a kink at T » T c (see, for 
example, Figs. 7 and 9), play only a minor role in the 
resistivity "maximum" at T c in optimally metallic GaM- 
nAs with most of the peak structure arising from the 
screening properties of ionized impurity scattering. The 
second generic experimental feature in Fig. 15, the slow 
decrease of p(T) for T > T c , cannot be explained at allby 
spin disorder scattering since spin disorder should remain 
large in the paramagnetic system (T > T c ) and certainly 
should not decrease with increasing T . Our Boltzmann 
theory provides a natural explanation (Figs. 13 and 14) 
for the decreasing p(T > T c ) as arising from the energy 
averaging of the relaxation time (Eq. 19), i.e., p(T) de- 
creases with increasing T simply because the holes move 
"faster" at higher temperatures (i.e. increasing kinetic 
energy with increasing T). This decreasing p(T) with in- 
creasing T (> T c ) also shows that our neglect of phonon 
scattering in the transport theory is a valid approxima- 
tion since phonon effects will always increase p(T) with 
increasing T. p(T) will increase again when phonon scat- 
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tcring starts dominating over ionized impurity scattering 
at much higher temperatures. The relative lack of impor- 
tance of phonon scattering in DMS systems arises from 
their very strong charged impurity resistive scattering ef- 
fects as reflected in very small sample mobilities. 

VII. CONCLUSION 

In this paper we investigate the transport properties 
of the diluted magnetic semiconductors using dynamical 
mean field theory and Boltzmann transport theory. We 
have shown that the resistivity depends strongly on the 
system parameters, i.e., exchange coupling, carrier den- 
sity, doping, and temperature. The resistivity drop with 
decreasing temperature in the ferromagnetic state can 
be partially explained by the screening theory for metal- 
lic samples. The parameter dependence of the resistiv- 
ity contains important information about the physics of 
diluted magnetic semiconductors. We find that in the 
strong exchange coupling regime the spin disorder scat- 
tering and the formation of the bound state in the impu- 
rity band compete to produce an unusual behavior in the 
temperature dependent resistivity. We also show that in 
the weak coupling regime the occupation of the minor- 



ity spin band is critical to the scattering mechanisms, 
and substantially reduces the resistivity because the re- 
pulsive interaction between local moments and "wrong- 
spin" carriers suppresses the carrier amplitude at the im- 
purity site, reducing the effective carrier-spin coupling. 
Our Boltzmann transport theory for charged impurity 
scattering is good qualitative agreement with the exist- 
ing DMS experimental data, showing that transport in 
DMS GaMnAs system may very well be dominated pri- 
marily by screened ionized impurity scattering (with spin 
disorder scattering playing only a minor secondary role), 
at least for the optimally doped metallic GaMnAs sam- 
ples. We have completely neglected detailed band struc- 
ture complications (e.g. spin-orbit coupling in the va- 
lence band) in our theory. These effects are certainly 
very important, but our interest in this paper is the de- 
velopment of a conceptually coherent qualitative theory 
for DMS transport identifying the main transport mech- 
anisms, and as such we have neglected all nonessential 
complications. 
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